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УДК 621.9.048.6 
Д.Н. МАЛЕЕВ, В.Х. АЛЬ-ТИББИ, Д.А. ЧИЛИКИН 


ОПТИМИЗАЦИЯ ПРОЦЕССА ЭЛЕКТРОАКУСТИЧЕСКОГО НАПЫЛЕНИЯ 
ПО КРИТЕРИЮ МИКРОТВЕРДОСТИ 


Исследуется влияние акустических параметров (амплитуда ультразвуковых колебаний А, частота ультразву- 
ковых колебаний Г и угол между осью волновода и направлением вектора продольно-крутильных колеба- 


НИЙ @ ) процесса электроакустического напыления защитных покрытий на качество покрытия (| микротвер- 


ДОСТ, Ъ). На базе проведенных исследований предложена математическая модель зависимости микротвердо- 
сти покрытия от акустических параметров процесса и проверена адекватность полученной модели. 
Ключевые слова: акустические параметры, амплитуда ультразвуковых колебаний, частота ультразвуковых 
колебаний, продольно-крутильные колебания, математическая модель, адекватность модели. 


Введение. Одним из перспективных методов нанесения токопроводящих износожаростойких за- 
щитных покрытий является метод электроакустического напыления - ЭЛАН [1,2]. Задача форми- 
рования качественного покрытия с высокими требованиями к шероховатости и микротвердости 
поверхностного слоя в процессе ЭЛАН является одной из важнейших при дальнейшем развитии 
этого метода для поверхностно-упрочняющей обработки деталей машин и формообразующего 
инструмента. Многовариативность выбора параметров процесса и их несомненное влияние на фи- 
зико-механические свойства покрытия, необходимость исследования физической картины процес- 
са требует проведения статистических исследований с целью обоснования режимов обработки и 
оптимизации показателей качества покрытия. Комплекс таких параметров существенно влияет на 
износостойкость и долговечность деталей машин и на другие их эксплуатационные характеристи- 
ки. В задачу исследований входило изучение влияния акустических параметров процесса ЭЛАН на 
качество наносимого покрытия. 

Методика исследований и некоторые априорные сведения. Исследование процесса ЭЛАН, 
безусловно, требует изучения влияния не только акустических, но и электрических параметров 
процесса на качество покрытия. Тем не менее, такая задача является слишком объемной в рамках 
одного исследования, поэтому в данной работе мы решили сосредоточиться для начала только на 
акустических параметрах процесса. Напряжение разрядного импульса как электрический пара- 
метр не претерпело изменений. По результатам предыдущих исследований [3,4] было выявлено, 
что удовлетворительную шероховатость покрытия (сопоставимую по классу с шероховатостью 
обрабатываемой поверхности) можно получить при работе на следующих усредненных по вели- 
чине режимах: напряжение разрядного импульса У=32 В (при напылении твердого сплава ВК6М 
на сталь 45), емкость разрядного конденсатора С=12000 мкФ, частота подачи импульсов 


Г. = 20 кГц. При этом скорость перемещения электрода составляла И =0,04 м/мин. 


Для разработки математической модели процесса был применен полнофакторный экспе- 
римент (ПФЭ), позволяющий получить математическую модель, выбрать контролируемые пара- 
метры и определить движение по градиенту, т.е. направление, при движении по которому пара- 
метр оптимизации уменьшается или увеличивается быстрее, чем в любом другом направлении. 
Были запланированы три вида статистических исследований, в каждом из которых фиксировалось 
три важнейшие характеристики качества покрытия: отношение микротвердости покрытия к мик- 


ротвердости подложки Нмп /Нмм, шероховатость покрытия К’ , мкм, толщина покрытия Н, мкм. 
В качестве параметра оптимизации у, , было взято отношение Н,п /Н,м, которое вычисля- 


лось после замеров на микротвердомере модели ПМТ-3, в соответствии с требованиями 
ГОСТ 9450-73. 
В качестве переменных (факторов) были выбраны следующие акустические параметры 


ЭЛАН (рисунок): х, - амплитуда ультразвуковых колебаний А (мкм), х, - частота ультразвуковых 
колебаний (кГц), х, - угол между осью волновода и направлением вектора продольно- 
крутильных колебаний сх (градусов). 
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Крутильная 
составляющая УЗК 


Продольная 
составляющая УЗК 


Схематическое изображение угла @ между осью волновода 
и направлением вектора продольно-крутильных колебаний 


Для описания процесса была выбрана линейная модель уравнения регрессии 
у=Б бл Вх лот лотах, (1) 
где у - параметр оптимизации; В, -— параметры модели; х,- факторы процесса ЭЛАН; { — номер 
параметра оптимизации, 1 =1,2,3.. 
Границы изменения факторов х,, х,, х, определены из установочных исследований [3,4] 
и соответствуют данным, приведенным в табл.1. 

















Таблица 1 
Значения исследуемых факторов 
Уровни Факторы 
А, мкм Х, кц @ ‚ град 
Верхний 17 22,8 82 
Нижний 4 16 10 
Основной 10,5 19,4 46 




















Интервал варьирования 6,5 3,4 36 
Планирование и проведение эксперимента. Был построен план матрицы планирования пол- 
нофакторного эксперимента типа 5. ‚ где 3 - число факторов. Кроме того, для обеспечения рас- 


четов ввели фиктивную переменную х, , которая принимает во всех опытах значения +1, а так- 
же графы взаимодействия факторов [5]. Данные сведены в табл.?2. 
































Таблица 2 
3 
Матрица планирования ПФЭ 2 
Номер точки Значения факторов Комбинации произведений факторов 
плана 
№ х Х› Хх хХ› хх: Хх: ХХХ 
1 + Е Е ь + + + - 
2 + + - - - - + + 
3 + - + р - + - + 
4 + + + - + - - - 
5 + ы ы + + ы - + 
6 + + Е + - + - - 
7 + - + + - - + - 
8 + + + + + + + + 
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Перед реализацией плана эксперимента опыты, предусмотренные матрицей, проводились 
в случайной последовательности. Порядок реализации опытов, выбранный по таблице равномер- 
но-распределенных случайных чисел, представлен в табл.3. 




















Таблица 3 
Рандомизация опытов 

Порядок Номер опыта 

реализации 1 2 3 4 5 6 7 8 
Серия 1 1 6 5 8 2 7 4 3 
Серия 2 4 2 7 5 3 1 6 8 
Серия 3 2 7 4 3 1 8 5 6 
Серия 4 8 1 6 7 5 2 4 3 
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Нужное значение факторов устанавливали поочередно в каждой точке плана матрицы в 
порядке реализации текущей серии эксперимента в соответствии с матрицей планирования и 


Н 
ТИ 
поддерживали эти значения постоянными в течение опыта. Параметр оптимизации У ( Е ) 
ИМ 





фиксировали по окончании опыта и заносили в табл.4. 
Таблица 4 





Е 
о Г! 
Результаты измерений параметра Ст ) 
































им 
Номер опыта Результаты эксперимента в сериях 

\11 \12 \13 \14 

1 1.79 1.91 1.73 1.86 

2 1.61 1.54 1.58 1.46 

3 1.77 1.81 1.79 1.84 

4 1.42 1.69 1.39 1.64 

5 1.58 1.41 1.60 1.42 

6 1.32 1.17 1.24 1.21 

7 1.49 1.27 1.41 1.35 

8 1.34 1.17 1.28 1.21 























Последующая обработка экспериментальных данных проводилась на ЭВМ, для чего была 
реализована программа в среде математического программирования Ма! Сад. 


Н ип 





Исследование зависимости = Г(А,у,@). Воспроизводимость измерений проверяли пу- 
ИМ 

тем подсчета среднего показателя параметра оптимизации процесса при реализации параллель- 

ных наблюдений, также была рассчитана выборочная дисперсия воспроизводимости параллель- 

ных наблюдений для каждой точки плана матрицы. Однородность дисперсий проверили по кри- 

терию Кохрена. 

Определили значимые коэффициенты регрессии и их значения с использованием крите- 
рия Стьюдента. Расчетные значения коэффициентов представлены в табл.5. В математическую 
модель процесса включили только значимые коэффициенты. Получили уравнение регрессии в 
виде 

У =1.50625 —0.11437. ХЛ - 0.16437 . ХЗ. (2) 
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Таблица 5 
Коэффициенты уравнения регрессии 
Коэффициент Значение { - критерий Значимость 
Ь, 1,50625 96,55 Да 
Ь, -0,11437 7,275 Да 
Ь -0,02062 1,28 Нет 
2 
Ь -0,16437 10,51 Да 
3 
Ь 0,02125 1,35 Нет 
12 
В 0,015 1 Нет 
Ь -0,00625 0,4 Нет 
23 
Ь 0,01312 0,83 Нет 
123 




















Проверка адекватности модели. Адекватность модели проверили по формуле Фишера. Таб- 
личное значение ГР =2,62 (при принятом количестве степеней свободы и уровне значимости). 
Поскольку Р’.(Е то гипотеза адекватности принимается. Были найдены коэффициенты 


расч табл! 


уравнения регрессии в действительных значениях. После простейших преобразований получили: 


Н,п 
—^ =1,901-0,0174-0,0046“ , 
им 
Так как в уравнении регрессии все коэффициенты, стоящие при х>, незначимые, то дан- 
ный фактор можно не учитывать. Фактор х, (амплитуда ультразвуковых колебаний А (мкм)) ока- 


зывает большее влияние на показатель параметра оптимизации по сравнению с фактором Х. 


(угол между осью волновода и направлением вектора продольно-крутильных колебаний сх (гра- 
дусов)). Это следует учитывать при контроле процесса. Следовательно, на микротвердость по- 
крытия наибольшее влияние оказывает амплитуда ультразвуковых колебаний. Существенен так- 
же параметр угла наклона вектора продольно-крутильных колебаний. Частота колебаний влияет 
незначительно. Знаки коэффициентов уравнения регрессии отрицательны, что говорит о том, что 
увеличение микротвердости можно достичь уменьшением параметров А и @. 

Методика оптимизации. Оптимизация проводилась методом крутого восхождения и начиналась 
от основных уравнений значимых факторов х,„, и х, „.Факторы изменялись пропорционально 


величинам коэффициентов регрессии с учетом их знаков. План крутого восхождения учтен при 
расчете в среде математического программирования Ма{1 Сад. 

С целью проверки полученной модели оценили показатель параметра оптимизации по по- 
лученному уравнению регрессии, используя рассчитанные значения факторов в плане крутого 
восхождения. Это необходимо не для предсказания выхода процесса, а для определения направ- 
ления движения. Данные "мысленных опытов " представлены в табл.6. 





























Таблица 6 
В Н 
"Мысленный эксперимент" для уравнения И =1,901-0,0174-0,0046а 
им 
Номер шага крутого восхождения А, мкм © «Мысленное значение» Е /Н и 

1 9.5 38.2 1.558 
2 8.6 31.4 1.605 
3 7.7 24.6 1.652 
4 6.8 17.2 1.702 
5 5.9 10 1.751 
6 5.07 2.8 1.8 

7 4.24 4.4 1.846 
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Полученный "выход" показывает, что направление движения выбрано верно - микротвер- 
дость покрытия возрастает. Вместе с тем показано, что нет необходимости в проведении опытов 
1-3, так как эффект варьирования параметров в указанных областях незначителен. 

Реализация оптимальных режимов для микротвердости покрытия. В соответствии с пла- 
ном крутого восхождения были проведены дополнительные эксперименты по вышеуказанной ме- 
тодике. Результаты экспериментов представлены в табл.7. Приведенные данные показывают, что 


экстремальное значение параметр Н „/Н м принимает в точке 5 при следующих оптимальных 


режимах: 4 =5,9 мкм, & = 10°. 




















Таблица 7 
Результаты дополнительных экспериментов 
Номер шага А, мкм © В ИН м 
3 7,7 24,6 1,91 
4 6,8 17,2 2,08 
5 5,9 10 2,15 
6 5,07 2,8 2,03 














Выводы. При постоянных параметрах напряжения разрядного импульса, емкости разрядного 
конденсатора и скорости перемещения электрода: 

1) на микротвердость покрытия наибольшее влияние оказывает амплитуда ультразвуко- 
вых колебаний; 

2) существенен также параметр угла наклона вектора продольно-крутильных колебаний; 

3) частота колебаний влияет незначительно; 

4) знаки коэффициентов уравнения регрессии отрицательны, т.е. увеличение микро- 


твердости можно достичь уменьшением параметров А и ©. 
5) оптимальные значения акустических параметров ЭЛАН по критерию твердости: ам- 
плитуда ультразвуковых колебаний А =12 мкм, частота ультразвуковых колебаний /=20 кГц, 


угол между осью волновода и направлением вектора продольно-крутильных колебаний сх =45°. 
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О.М. МАГЕУЕМ, \М/.Н. АТ-ТТВВТ, О.А. СНТЕТКТМ 


ОРТТМТРАТТОМ ОЕ ТНЕ ЕТЕСТКО-АСОЦ$ТТС 5РИТТЕВТМС РКОСЕ$З$ 
ВУ ТНЕ МТСКОНАКОМЕ$$ СВЕТТЕВТОМ 


ГпЯчепсе оГ асои$Ис рагатеег$ (атрЖиае оЁ иЁгазопс озсИаНоп$ А (тст), Гедиепсу оЁ иЁгазотюс 
о5сИаНоп$ (КП7), апа а согпег Бебмееп {Пе ах! о! {Пе ммауе дуае апа {Пе уесог атесйоп ог оподКч- 
та!-//$Чпта озсИаНоп$ (4едгеез)) ргосез5 оЁ ап ее го-асоц$Ис дизИпа оЁ {Пе рго{есфуе соа!пд$ оп 
дчаШу ог соаНпа (пусго Вагапез$) 15 пуезНдае4. Оп {Пе Ба$5 оЁ {Пе ехесше4 гезеагсй {Пе тае- 
таЧса! подЕ! оЁ дерепдепсе ог пусго Пагапез$ оЁ а соа!пд гот асои$Ис рагатецег$ оЁ ргосез$ 15 де- 
уеюре4 апа адедчасу оё {Пе гесеме4 тоде| 15 сПескеч. 

Кеу\мгог4$: асоц$Ис рагатеег$, атрийе оЁ иЁгазопю озсИа{оп$, Недиепсу о! ийгазотс озсИабопз, 
опдкипа!-бл$Нпоа озсИаНоп$, па ета{са! по@Е|, адедиасу о {Пе тодЕ|. 
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